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Abstract 

Parametric  spline  corvee  an  d  smfacas  a re  typically 
constructed  so  that  some  number  of  derivatives  match 
where  the  curve  segments  or  surface  patches  abut.  If 
derivatives  op  to  order  n  are  continuous,  the  segments 
or  patches  are  said  to  meet  with  C",  or  order  pa¬ 
rametric  continuity.  It  baa  bees  shown  previously  that 
parametric  continuity  Is  soft  dent,  bat  not  nscsssary,  for 
geometric  smooth  new. 

The  geometric  measures  of  sad  toageot  and  cur- 
eetarv  vectors  for  corvee,  and  tangent  plane  and  Dnpin 
indicetrix  for  surfaces,  have  been  seed  to  dsflae  first  and 
aacond  order  geometric  ceatiaaitg.  In  this  paper,  we 
extend  tbe  notion  of  geometric  continuity  to  arbitrary 
order  a  fC*)  for  corvee  and  surfaces,  and  present  an 
Intuitive  development  of  constraint  equations  that  are 
accessary  and  sufficient  for  it.  The  constraints  fin  own 
as  the  Bets  constraints)  result  from  n  direct  application 
of  the  univariate  chain  rule  for  come  and  the  bivariate 
chain  rule  for  surfaces,  lb r  Cist  sad  aacond  order  con* 
tinuily,  the  Beta  constraints  an  equivalent  to  requiring 
continuity  of  the  geometric  measures  described  above. 

The  Bets  constraints  provide  for  the  latrodnctioa 
at  quantities  haown  as  dope  parameter*.  If  two  cone 
segments  an  to  meet  with  C“  continuity,  n  shape  pa- 
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mnstars  map  be  Introduced,  fbr  surfaces,  the  ass  of 
the  constraints  for  G"  continuity  provides  for  the  intro¬ 
duction  of  n(n  +  S)  shape  functions,  defined  along  the 
boundary  between  two  surface  patches,  fbr  polynomial 
splines,  the  ess  of  the  Beta  constraints  allows  greater 
flexibility  through  tbe  shape  parameters  without  rais¬ 
ing  tbe  polynomial  degree. 

The  approach  we  tabs  Is  Important  for  asveral 
reasons:  first,  it  generalises  geometric  continuity  to 
arbitrary  order  for  both  carves  and  surface*.  Second, 
it  shows  the  fundamental  connection  between  geometric 
continuity  of  curves  and  that  of  suifocus.  Third,  dus  to 
tbs  chain  nils  derivation,  constraints  of  any  order  can 
be  determined  more  easily  than  using  derivations  based 
exclusively  on  geometric  measures. 

Riaumt 

Lee  combes  et  surfaces  paramdfriquee  I  base  ds 
splines  son!  gdodrafemenl  coastruites  de  fog***  icequ’un 
certain  a  ombre  de  ddrfvd as  coincident  sux  raccorde- 
ments  antra  las  tree  de  combe  on  las  carresux  de  sur- 
foce.  Lonqu’addiibnaelhment  Ise  a  premMrm  ddrivdm 
noat  continues,  ke  arcs  on  las  carreaux  as  rencontreot 
a  vac  eontinuiU  pemmitrigae  C",  oe  d’ordra  n.  0  a  ddji 
dtd  ttabli  qus  la  coatinuitd  paraodtrique  set  suffisants  h 
i’obtcnfion  d*ua  image  gdoi mdtrique,  mais  qu'elie  n’est 
pas  adcassaire. 

Las  premier  et  deuxMme  ordree  de  continuity 
gdomdtriqu*  sent  gdadraiemenl  ddflais  h  i’aide  de 
meeurae  gdomdtriquas  tela  Is  ssetser  Mage  at  uatfatri 
et  Is  eeetcer  ds  couriers  dans  Is  cas  das  combes,  aiasi 
qus  Is  pie »  leaf*  at  et  ISndiceirix  de  Depin  dans  h  cat 
dee  surfaces.  Dans  cat  art  ids,  nous  gdadralisoas  la  to- 


Graphic*  Interface  ’85 


tfca  de  eoatlauM  gdoaadtriqoe  A  m Import*  qool  or dn  m 
(C*J  aooeJ  bha  pour  ho  o ourbm  quo  poor  ha  surface*. 
Horn  pedoentoae  dgakmmt  aa  ddrehppmamt  btuHif 
dm  dquatioae  do  eoatraiat*  adeomalrm  ot  mihirtw 
Cm  e oatniatee,  go*  too*  appeioae  Im  emfmtafes-lete, 
ddeoohat  diroetemeot  dm  rkgJm  do  cfcalbe  A  im  vuri- 
obk  poor  ho  eourbm  ot  A  d*ux  nriablm  pour  ho  oar 
to cm.  Poor  ho  prmohr  ot  oacoad  ordrm  do  coatbuitd, 
Im  eootnlatm-b* la  soot  dqahrahatm  A  h  coatbuitd  dm 
am tm  gdomdtriqum  dderltm  d-dmmto. 

Im  coatrafatee-bets  otrmt  I’occaekn  dlatrodolro 
eertaiam  quaatitde  eooaum  oooo  h  oom  do  paramitrea 
do  form* i .  8i  doux  arm  do  eoorbo  dohomt  m  taceotdm 
iwc  eoatiauMG",  m  parambtrm  do  formmpmmmt  Ha 
latmduito.  Poor  hooortoeoo,  I'otllioatha  dm  CMtnbto 
d*  eoatlaoM  C*  pormot  dlatrodolro  «(•+*)  function* 
do  form**,  hoqoolho  ooot  ddiaho  h  hog  dm  Hmltm 
eommuammtrthooarAcmdodooacanoouamheyoou 
Doao  h  cm  dm  apliam  pofyacenlaax,  I’otibootioo  dm 
ooatntotm-boto  pormot  dm  InttiVM  accrue  pic*  on 
parambtrm  d*  hrmm,  oooo  poor  Mtul  uimittr  h 
dagrd  do  pofyo&m*. 

Colt*  opprock*  mt  Important*  poor  molatm 
raisons.  Promhromoat,oUogda4nBoohaohoadoeooti- 
aoM  gdorndtriqum  iu  ordrm  qaokoaqam,  lltut  poor 
ho  eourbm  quo  pour  ho  oarhcoo.  Urndto*— I,  0O0 
amt  oa  dvideac*  la  aimllaritd  toadamoatah  mt n  h  com- 
tlaultd  gdomdtriquo  dm  eourbm  ot  ealh  dm  oortaeaa.  FU 
aahmont,  ha  righo  do  chalao  hcllltmt  la  ddtermbatioa 
do  coatraiaim  d’ordro  qookoaqao,  eomparnthommt  A 
eoqu’aurahat  ponahdm  ddrfratioaa  booh*  oaiqaoaamt 
our  dm  maouim  giomdtriqaoa. 

KEYWORDS:  pomrtric  moAtlliaf,  coatbaity,  pma- 
tirk  cimi,  ptniMtrk  ourlacm ,  shape  p nrnmoton. 

1.  introduction 

Cwwd  mb  defined  or  youofutedby  parumitriaohom 
[ourlacm  will  b«  addressed  b  lection  I).  A  oaooori- 
oto  (om  variable)  parametrisation  b  a  fiactlo*  such  M 
«(•)  •  (Jr(w),F(o)),  where  tbo  domain  parameter  ,  k 
•flowed  to  rang*  om  ootn*  btanl  ]«•,  «i  J.  hr  a  given 
value  of  «,  tbe  function  ({«)  cm  bo  thought  of  m  local- 
b<  »  partkk  b  Euclidean  two-spae*.  Aa  «  b  Increased 
om  tbo  bteml,  tbo  partick  tnww  a  path  defined 
by  q,  t raciaf  out  a  com  b  tbo  procoM  (ooo  Figure  1). 
If  l«o,  Ml]  b  tboofbt  of  M  aa  oriented  Hat  augment,  then 


4  eaa  bo  tboafbt  of  m  a  deformation  producing  as  ori¬ 
ented  eon*.  Tbo  Ant  dorintni  rector  q(l>  repnoento 
tbo  velocity  of  tbo  partkk  (b  fonoral,  «o  denote  tbo  i*k 
derivatiro  of  a  univariate  foactioo  by  onperociipt  (»))• 
Tbo  rolodty  b  a  rector  quantity  and,  aa  ouch,  contains 
Information  about  orientation  and  rate,  or  epood.  Tbo 
aocond  derivatiro  rector  q(,)  roproooata  tbe  acceleration 
ot  tbe  partick,  oo  it  too  contains  information  about  tbo 
(change  of)  rate.  Thus,  a  parametrisation  contains  in- 
fonnatioo  about  tbo  geometry  (tbo  obape  or  image  of 
tbe  cum),  tbe  orientation,  nnd  tbo  rate. 


Ilgen  1.  Tbo  nnioariatt  parometriaation  q  generate* 
on  oriented  enroe  bp  deformation  of  the  oriented  line 
oofmsMf  |«o,«i). 


Figure  2  show  tbe  curves  generated  by  three  differ¬ 
ent  parametrisation*.  Tbe  obape  of  tbe  curreo  b  identi¬ 
cal;  they  differ  only  b  orientation  and  rate.  Curreo  (a) 
and  (b)  bar*  tbe  earn  orientation  at  each  point,  but  the 
rateo  differ.  Tbe  cum  labelled  (c)  differ*  from  (a)  and 
(b)  b  orientation  and  rate.  If  a  cum  b  defined  to  be 
■imply  tbo  geometry  property  of  a  parametrisation,  one 
vould  conclude  that  figures  (a),  (b),  and  (c)  represent 
ofwoolcnf  curves.  Wo  will  refer  to  this  as  the  C  model 
el  a  cum.  Aaothsr  possibility  b  to  consider  tbo  geom¬ 
etry  and  orientation,  which  wo  wiD  call  tbo  GO  model 
Using  tbo  OO  model,  ooo  would  say  that  (a)  and  (b) 
are  equivalent,  but  (c)  b  different.  The  last  poMibil- 
ffey  wo  will  consider  b  tbo  GOR  model,  where  geometry, 
orientation,  and  rate  me  nil  relevant  to  the  definition 
ot  a  curve.  Using  this  model,  uo  pair  of  the  curves  b 
equivalent. 

In  recent  yean,  heavy  nee  has  bees  made  of  piece¬ 
wise  parametric  functions  known  aa  parametric  apbin e*. 
Spline  came  ere  typically  constructed  by  stitching  to¬ 
gether  nahrariato  parametric  functions,  requiring  that 
some  number  of  derivatives  match  at  each  joint  (the 
points  where  tbe  cum  segments  meet).  If  n  derivatives 
agree  at  a  given  jobt,  tbe  parametrieatione  there  ere 
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Figure  I.  Back  of  ike  nmi  okoot  hot  tkt  tom 
image;  they  onlg  differ  in  trinMi'M  mod  rati.  Orien¬ 
tation  it  indicated  If  arrowheads  to i  rate  it  radfrsted 
If  outon  tangent  to  tkt  tones. _ 


•tud  to  nut  with  »tt  order  pa rtmttrie  ceatinstty  (C* 
contbuity  for  short). 

Wc  maintain  that  tbs  choies  of  a  particular  nodd 
for  a  curve,  and  hence  the  choice  of  how  the  cam  seg¬ 
ment*  are  stitched  together,  eh  on  id  be  application  de¬ 
pendent.  Fbr  instance,  if  a  spline  is  being  need  to  define 
the  motion  of  an  object  in  an  animation  system,  the 
GOR  model  is  most  appropriate  sine*  the  orientation 
and  rate  are  of  importance,  la  this  type  of  applica¬ 
tion,  parametric  continuity  is  required  to  maintain  the 
smooth  nee*  of  the  rat*  properties.  In  other  words,  par 
ramctric  continuity  will  ensure  that  the  object  will  moot 
smoothly. 

However,  in  CAGD  the  rate  aspect  of  a  parame- 
trisatioo  is  often  unimportant.  Coarider  for  example 
the  use  of  aplines  to  describe  numerically -con  trolled  cut¬ 
ters.  It  may  be  necessary  to  specify  uniquely  th*  direc¬ 
tion  of  the  cutter  at  each  point  on  th*  path,  but  th* 
speed  of  th*  cutter  may  depend  upon  th*  hardness  of 
the  material  being  cut.  For  this  type  of  application,  the 
GO  model  is  most  suitable,  but  parametric  continuity  is 
overly  restrictive  tine*  it  places  emphasis  on  irrelevant 
rat*  information.  Many  other  applications  b  CAGD 
require  only  th*  G  model,  but  it  seams  difficult  to  de¬ 
velop  a  useful  formalism  without  th*  structure  provided 
by  orientation.  We  will  therefor*  adopt  th*  GO  model, 
and  develop  an  appropriate  measure  of  continuity,  one 
based  based  only  on  th*  geometry  and  orientation  prop¬ 
erties;  we  refer  to  this  as  geometric  ceatinudf 

It  has  recently  coma  to  our  attention  that  many 
author*  have  Independently  defined  this  kind  of  con¬ 
tinuity  of  first  and  second  order  (which  we  denote  by 
C‘  and  CP,  respectively)  for  curves  and/or  surfaces  us¬ 


ing  geometric  means.  Fbr  curves,  Fbwler  It  Wilson1*, 
8abin,T,  Manning1*,  Phan  It  Pratt*,  and  Baraky*  each 
independently  defined  first  order  continuity  by  requir¬ 
ing  that  the  salt  tangent  oacton  agree  at  th*  jobts.  To 
achieve  second  order  continuity,  both  the  unit  tangent 
and  curvature  vectors  were  required  to  match.  Niel¬ 
son's  v-eplbc14  possesses  a  similar  Mad  of  contbuity. 
Those  geometric  measures  essentially  ignore  th*  rat*  in¬ 
formation  by  Normalising  the  parametriMtioa  before 
determining  smoothness. 

Fbr  surfaces,  H  is  common  to  require  matching  of 
fsafcuf  planet  tor  first  order  geometric  contbuity  (cf. 
fiabb1*  and  Veroa  «t  aP°).  Ibr  surfaces  of  second  or¬ 
der  geometric  contbuity,  Veroa  at  al  and  Kahmann1* 
require  contbuity  of  normal  csrswfsre  b  every  direc¬ 
tion,  at  every  point  on  th*  boundary  shared  by  the  con¬ 
stituent  surface  patches.  As  Vbron  et  al  and  Kahmann 
each  show,  this  is  equivalent  to  requiring  that  th*  Dopio 
iadusMt  (cf.  DoCarmoT)  of  each  patch  agree  at  the 
boundary  curve.  The  Dupb  bdicatrix  is  a  measure  of 
curvature,  but  the  curvature  properties  of  surfaces  are 
sufficiently  complex  that  they  cannot  be  characterised 
by  something  as  simple  as  a  scalar  or  a  vector. 

Although  the  geometric  approaches  described 
above  are  convenient  and  btuitiv*  for  first  and  second 
order  contbuity,  n  more  algebraic  development  is  bet¬ 
ter  salted  for  th*  extension  to  contbuity  of  higher  order. 
The  approach  we  take  is  based  on  rapnrametriMtioa  — 
the  process  of  obtaining  a  new  panmetrbatioa  given  an 
old  on*.  Undo-  the  GO  model,  rsparnmetrisation  may 
change  rate,  but  not  geometry  or  orientation.  By  allow¬ 
ing  reparametrbntbn  before  making  n  determbation  of 
continuity,  the  rate  aspect*  of  p  ranetriiations  may  be 
ignored.  Alternately  stated,  j  approach  is  baaed  on 
the  following  simple  idea: 

Pi;  Don’t  base  contbuity  on  th*  parametrisations  at 
hand,-  reparamstrise,  If  >ac senary,  to  obtain  param- 
otriaatloaa  that  meet  with  parametric  continuity. 
If  this  caa  bs  doom,  th*  original  pararaetruathaa 
moot  alto  asset  amootkfy,  at  least  b  a  geometric 


The  above  concept  is  not  n  new  one;  similar  prin¬ 
ciples  have  been  dieenmed  by  fhria*  and  Venn  et  al**. 
What  is  sew  is  the  ass  of  th*  principle  to  construct  con- 
strabt  equations  (known  as  the  Beta  constraints)  that 
are  necessary  and  sufficient  for  geometric  contbuity  et 
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arbitrary  order  for  both  cam*  lid  —hm,  * 

The  Beta  constraint*  generalise  the  penmetric  con¬ 
tinuity  constraints  through  the  Introduction  of  freely 
variable  quantltiee  called  shape  parameter*.  Once  the 
Beta  constraint*  are  determined  for  a  given  order  of  con¬ 
tinuity,  they  may  be  need  in  place  of  the  parametric  con¬ 
tinuity  constraints  when  building  ^liaes,  thereby  ob¬ 
taining  increased  flexibility,  Fbr  instance,  if  the  C4  con¬ 
straint*  are  replaced  with  the  C*  constraints  in  the  uni¬ 
form  cubic  Bopline",  the  cubic  Beta-spline  results1 -1. 
The  cubic  Beta-spline  is  an  appnwimsfinf  spline  tech¬ 
nique  that  possesses  two  shape  parameters;  an  caUr- 
poloting  technique  is  described  in  DeBoee  At  Bnraky*. 
Faux  Ac  Prntt*  and  Fbria*  nee  the  extra  freedom  allowed 
by  geometric  continuity  to  pines  Bit ier  control  Metises. 

As  important  aspect  of  thane  techniques  b  that  the 
additional  flexibility  of  geometric  continuity  is  added 
without  increasing  the  degree  of  the  polynomials.  Thin 
is  particularly  important  for  algorithms  that  manipulate 
the  spline.  Fbr  instance,  the  complexity  of  Sederbergb 
algorithm"  to  intersect  two  polynomial  carves  of  degree 
d  grows  at  least  as  Cut  as  d*.  Substantial  savings 
can  therefore  be  bad  by  mlnlmlaing  the  degree  of  the 
polynomials  involved. 

la  the  remainder  of  this  paper,  we  extend  the 
notion  of  geometric  continuity  to  arbitrary  order  n  (G*) 
and  show  (in  n  nonrigorous  way)  that  the  derivation 
of  the  Beta  constraints  results  from  n  straightforward 
nee  of  the  univariate  chain  rale  for  carves  and  the 
bivariate  (two  variable)  chain  rale  for  surfacee.  fbr 
n  more  complete  treatment,  the  reader  b  referred  to 
Banky  Ac  DeRoee*  and  DeRoee*. 

3.  Geometric  Continuity  for  Cnrren 

We  begin  the  study  of  geometric  continuity  for 
curves  by  examining  the  rsparametrisatioa  proem*. 
Two  parametrisatioBS  are  raid  to  be  00  eyw'enfraf  If 
they  have  the  same  geometry  and  orientation  is  the 
neighborhood  of  each  point.  Given  n  panmetriention 
q,  all  GO- equivalent  panmetrisatione  may  be  obtained 
by  function*!  composition.  More  specifically,  If  q(u) 
and  5(5)  are  GO-equhralent,  thm  they  an  related  by 
5(5)  -  q(e(5)),  for  some  appropriately  choeea  change 


of  pommeter  e(S)  (see  Figure  fl).  Since  q  and  5  must 
have  the  same  orientation,  •  must  be  an  increasing  func¬ 
tion  of  5,  implying  that  •  must  satisfy  the  orientation 
press  mag  condition  «<>>  >  a  Intuitively,  u(5)  deforms 
the  interval  (So, Si)  into  the  interval  [uo,«i]  without  re¬ 
versing  the  orientation  of  the  segment  fc,5||.  This  in 
tan  implies  that  q  and  5  will  have  the  same  geometry 
and  orientation,  bat  they  may  differ  In  rate. 


“0  «1 


Figure  S.  The  efusvefcaf  pammetriaaiione  q  sad  5 
arc  nlntid  fy  the  change  of  pommeter  w(5). 

A  oahrariate  parametrisatioa  is  regular  if  the  first 
derivative  vector  dose  not  vanish.  It  is  well  known  from 
differential  geometry7  that  regularity  is,  in  general,  es¬ 
sential  for  the  smoothness  of  the  moulting  curve.  We 
will  therefore  restrict  the  discussion  to  regular  parume- 
trisationa.  We  sow  ghee  n  more  precise  definition  of  G* 
continuity: 

Definition  It  6  (<o.ti)  and  q(u),u  €  |«o,s,| 

ho  two  parametriaatioae  tack  that  r((i)  ■  q(«o)  (nee 
figure  4).  These  parametriaatioae  meet  with  C"  cooti- 
aaJty  at  J  V  and  only  if  there  exist  GO-equ/vaJent  pa- 
rametriaatioae  f (?)  and  5(5)  that  most  with  C*  conti¬ 
nuity. 

Definition  1  is  simply  n  restatement  of  principle  Pi, 
but  to  practice  one  cannot  examine  all  GO-equhralent 
parametrisatioa*  to  an  effort  to  find  two  that  meet  with 
parametric  ceottoaity.  However,  H  is  possible  to  Cad 
conditions  on  r  and  q  that  an  necessary  and  sufficient 
for  the  emstence  of  GO-equhralent  panmetrisatione  that 
meet  with  parametric  continuity. 

Although  Definition  1  suggest*  that  both  r  and  q 
need  to  be  reparametrised,  it  is  possible  to  show  that 


Graphics  Interface  ’88 


r(<,)=q(u0)=y 


Flpn  4.  The  parametrisations  r(<)  on 4  q(«)  mmt  at 
the  common  point  J. _ 

Definition  1  holds  if  sad  only  if  them  exists  a  $  that 
meets  »  with  parametric  continuity.  la  other  words,  only 
oae  of  the  parametrisations  seeds  to  be  reparamstriaed 
to  determine  smoothness. 

We  will  ultimately  be  interested  ia  the  derivative 
properties  of  q.  The  viueenefe  chain  rale  allows  as  to 
express  derivatives  of  q  in  terms  of  the  derivatives  of  q 
and  *.  For  example,  the  first  derivative  ia  given  by 

q  a  45 

m*ti a  (*») 

tie  dm 

In  general,  the  Ith  derivative  of  q  can  be  srrittea  as  some 
function,  call  it  CJti,  of  the  first  i  derivatives  of  •  and 
q.  That  is, 

•">. -,•<*>).  {  ) 

We  are  actually  interested  is  qW  evaluated  at  its  le/t 
parametric  endpoint  5<>.  Thus,  derivathrss  of  q  and  ■ 
must  also  be  evaluated  at  their  left  endpoints: 

q(,‘(5o)  -  C*<(q“»(ue),-.,qW(uo), 

Since  «  is  a  scalar  function,  evaluating  one  of  its 
derivatives  results  ia  a  real  number.  In  particular,  1st 
sw(ii)  ■  Pj,  i  ■»  1, ...,«.  Equation  (14)  then  becomes 

^(Ss)  -  Cl<(q,,,(ue),  --,q<<,(«e), 

The  orientation  preserving  quality  of  n  implies  that 
Pi>0 


We  are  now  in  a  position  to  state  the  primary  mult 
of  geometric  continuity  for  curves,  tec  all  that  r  and  q 
meet  with  G*  continuity  if  q  can  be  reparametrised  to  q 
so  that  derivatives  of  w  and  )  agree.  That  is,  we  require 
that 

»*°(*a)  *  q<0(!o).  f  ■  !,—«».  (16) 

Positional  continuity  is  implicitly  assumed  (see  Fig¬ 
ure  4).  Substituting  equation  (14)  into  (II)  yields 

e,0(<i)  •  Cl<(q(,,(«e),  •  •  •  ,q,0(ne),  . 

Pi.-.P.)  ,“'*- 

CM) 

The  constraints  reeultiag  from  equation  (14)  are  the 
saiserisfe  Beta  constraint*  and  the  numbers  A , 
are  the  shape  parameters.  The  above  discussion  is  not 
n  proof  that  the  Beta  constraints  are  necessary  and 
sufficient  conditions  for  geometric  continuity,  but  such 
a  proof  can  be  constructed*''.  Thus,  if  equations  (16) 
are  satisfied  for  any  choice  of  the  fia,  subject  to  P,  >  0, 
then  the  coincident  curve  segments  will  meet  with  G“ 
continuity,  fbr  instance,  the  Beta  constraints  for  G* 
continuity  between  t  and  q  are 

r(,)(f|)«  Aq(i>(«o) 

»‘,|(f.)-«V,|W+Aq"lW 

rw(l.)  -  PI  qW(«o)  +  IPX Ps  q,,»(ne) 

+Aq,1,W  [ 

»w,(f.)  -  Pi q,4,(«e)  +  A  «W(«e) 

+  (4 PiPs  +  fi^)q(>,(«o)  ♦  Aq0,W' 

Although  equations  (17)  were  derived  using  the 
chain  rule,  the  first  two  are  identical  to  the  constraints 
resulting  from  a  geometric  derivation  of  unit  tangent 
and  curvature  vector  continuity'  ".  Thus,  our  approach 
reduces  to  previous  definitions  of  C1  and  C  continuity 
for  curves.  It  can  also  be  shown  that  Beta  constraints 
for  »**  order  continuity  are  equivalent  to  requiring 
continuity  of  the  first  n  derivatives  with  respect  to  arc 
length  *'*. 

When  constructing  a  spline  technique,  if  the  Beta 
constraints  are  used  in  place  of  the  parametric  conti¬ 
nuity  constraints,  new  freedom  is  introduced  through 
the  shape  parameters.  These  parameters  may  be  made 
available  to  a  designer  in  a  CAGD  environment  to 
change  the  shape  of  the  target  curve,  as  the  following 
example  shows. 


Buapl*  S.li  1b  dtBoutnta  Ikmcflkt  B«U 
constraints,  we  «IU  sketch  Ik*  construction  of  Ik*  go- 
oowtric  continuous  analogue  of  Ik*  uniform  quartk  B- 
spline  called  (naturally  enough)  Ik*  quartk  Bsta-spline. 

The/**  segment  of  tk*  qoirtic  B-epUne  i*  generated 

«,(•)-  £  ▼#***(•).  •  €  jO,  1]  (2J) 

*■-* 

where  Ik*  hsm  functions  £*(«)  or*  quartk  polynomials 
that  satisfy 

*&(»)- 4t°( 0). 

<-0,  (*A) 

1  se  "■  2,  ihj  I  ■ 

Tk*  **qn«Bo*  of  control  oorHeoi  Vy+*  comprise  a  control 
pofyfon. 

Biac*  tk*  dtrintin  properties  of  Ik*  baaia  functions 
are  inherited  by  q, ,  aquation  (2.9)  implies  that  tk*  cam 
segments  meet  with  C*  continuity.  Tk*  quartk  Beta- 
spline  i*  constructed  by  building  quartk  pofynomiak 
M«)  that  satisfy  tko  O*  constraints  instead  of  tk*  C* 
coast raiats  of  equation  (2.9).  That  is, 

-«<(»"’(•).  -.•C’w.  .  . ... 

‘  ‘U'U' 

(1.10) 

Eqnatioa  (2.10)  implies  that  Ik*  basis  functions  an 
dependent  upon  lb*  shape  parameter  values.  Changing 
a  shape  parameter  therefore  chaagas  the  shape  of  the 
resulting  curve  (see  Figure  S).  o 

S.  Geometric  Continuity  for  Surfaces 

la  this  Mctioa,  we  extend  the  notions  of  geometric 
continuity  to  surfaces.  Since  care  vas  taken  ia  Section  2 
aot  to  base  the  development  of  geometric  continuity  on 
concepts  (such  as  arc  length)  that  don’t  apply  to  sur¬ 
faces,  Ike  machinery  developed  for  anlvariat*  parame- 
trisatioas  can  be  readily  extended  to  bivariat*  parame- 
trteatioas. 

A  onrfoct  patch  k  deftaed  by  a  bivariat*  fuactioo 
each  a*  C(e,e)  -  (X(u1o),Y(n, vj.lfe,*)),  where  • 
aad  e  are  allowed  to  range  over  some  region  D  of  the 
•e  plane  (a*e  Figure  9).  Loosely  speaking,  a  onrfoct  k 
a  collection  of  surface  patches.  We  use  the  notation 
C,,-fl(e,v)  to  dcaote  the  i1k  partial  derivative  with 


Figure  B.  The  cerecs  sieve  share  the  same  control 
polygon,  and  all  h*se  fa  ■  1  and  fa  m  0;  they  differ 
only  in  the  sslne  of  fa.  The  top  enree  Uo  fa  *0,  Uu 
middle  cores  has  fa  «  20,  and  the  bottom  enree  has 

ft  -  100- _ 

respect  to  a,  aad  the  j*k  partial  with  respect  to  «.  Ia 
general,  a  superscript  (»,/)  denotes  the  i*k  partial  with 
respect  to  the  first  variable,  aad  the  j*  partial  with 
respect  to  the  second.  A  bivariate  parametrisatioa  each 
as  G  k  regular  if  the  first  order  partial*  (G*10*  aad 
Qi*-1))  are  liaaarly  independent;  we  will  deal  exclusively 
with  regular  parametrkationa. 

Ia  Section  1,  we  aaw  that  univariate  parametrka- 
tlooe  contain  iaformatioa  about  geometry,  orientation, 
aad  rate.  The  tame  k  true  of  bivariate  parametrkatioas. 
Orientation  can  be  defined  by  treating  D  as  an  oriented 
plane  haviag  a  "top  side*  and  a  "bottom  ride.*  G  can 
than  be  thought  of  as  deforming  the  oriented  plans  to 
produce  an  oriented,  or  two-ekled,  surface  patch.  The 
rate  information  enters  through  the  partial  derivatives 
of  the  parametrisatioa.  We  can  therefore  speak  of  the 
G,  GO,  aad  GOR  models  of  eurfacae.  Just  as  for  curves, 
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the  of  a  particular  modal  should  ba  application  de- 
paodaot.  Wa  will  adopt  the  GO  modal  for  two  naaona: 
tnt,  orientation  ia  necessary  is  applicationa,  each  aa 
rendering,  where  tba  two-aid adneaa  of  surfaces  ia  impor¬ 
tant,  aad  aacoad,  it  aaama  difficult  to  develop  a  aoafal 
lonnaliam  without  tba  atructun  provided  by  orienta- 
tioo,  especially  when  surfaces  am  allowed  to  btemect 
tbamsehras. 

We  aow  axamiaa  tba  reparametriaation  process  for 
surface  patches.  Two  bivariate  parametriaationa  an 
GO- equivalent  if  they  have  the  same  geometry  aad 
orieatation  ia  the  neighborhood  of  each  point  oa  tba 
surface  patch.  If  C(«,  •)  aad  C(«,  v)  are  GO-equivalent, 
then  they  an  related  by 


Just  aa  for  curves,  parametric  continuity  ia  appro¬ 
priate  for  tba  GOR  model  of  a  surface,  but  it  ia  not 
suitable  for  use  with  the  GO  model  dace  it  places  em¬ 
phasis  oa  irrelevant  rate  information.  The  determina¬ 
tion  of  continuity  can  be  made  insensitive  to  rate  by 
allowing  reparametriaation.  Thui,  wa  aay  that  P  and  G 
meet  with  C"  continuity  if  mud  only  if  then  exist  GO- 
equivalent  parametriaationa  P  and  C  that  meet  with  C* 
continuity. 


<t6 


C(*,e)-C(u(u,5),e(ul5)) 


where  the  functions  u  aad  e  aatiafy  the  orientation 
preserving  condition  ** 


u(tT,tT) 

v(u,v) 


w(I.S)v(0.t)  _  m(0.l)w(l,0)  >  0 


Ws  aow  examine  how  surface  patches  an  stitched 
together  with  parametric  continuity.  Referring  to  Fig- 
are  7,  P(e,l)  and  C(u,e)  meet  with  »**  order  param¬ 
etric  continuity  if  and  only  if  all  like  partial  derivatives 
of  order  up  to  a  agree  for  each  point  of  the  boundary 
carve.  That  is, 


Pigure  •.  G  end  G  era  CO-qeiWiel  parametrise- 
tiom  raiefad  fy  the  c  he  aye  of  poromttrimhon  dtttr- 
mitud  fy  «(«,  5)  end  v(5, v). 


P^yJ-C^y),  <  +  ,*- l.^n, 


when  evaluation  at  y  la  to  be  interpreted  aa  evaluation 
at  all  points  P  of  y. 


In  complete  analogy  with  curves,  only  one  of  the 
panmetrisations  actually  Beads  to  be  re  parametrised, 
implying  that  P  aad  G  meet  with  G*  continuity  if  and 
only  if  then  exists  a  C  such  that 


**  liedrte  fcwllier  eith  — UKerises  eatrehi 
•qaeStee  (».»)  m  the  Utah Urn  at  I he  cheap  at 


!“•%)«  eW(y), 


(<f.  D*Carmar) 


Once  again,  b  complete  analogy  with  curves,  the  bivari¬ 
ate  chain  rule  can  be  used  to  express  derivatives  of  G 
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•**» 


1b  terms  of  G.  1b  general,  the  i,  j**  partial  of  6  caa  ba 
expraeeed  a*  now  function,  call  H  at  the  partial* 
of  G, «,  ud  «,  ap  to  order »+/.  Statad  mathematically, 

G,,jl  -  CZij{GlkJ),»{kJ),t{kJ)),  (SJ) 

where  the  ifidicee  (4,1)  an  to  take  oa  all  positive  values 
each  that  l  +  Ui-f/. 

We  may  now  obtain  the  Bivariate  Bata  constraint* 
by  evaluating  (S.S)  along  the  boundary  carve,  followed 
by  eubotitntion  into  equation  (1.4)  to  get 

T^h)  -  C^(G<*  ')(y), 

tor  k  +  lmi  +  j  and  •  +  ;•  1,^*. 

Hie  equation*  leeulting  from  (S.6)  are  the  bivariate 
Beta  constraint*,  and  the  acalar  function*  and 

v(*.>)(^)  are  the  ahape  functions.  A  simple  counting 
argument  shows  that  it  is  possible  to  introduce  a(n + 8) 
shape  functions  when  two  patches  are  stitched  together 
with  G*  continuity. 

Just  as  the  univariate  Beta  constraints  can  sup¬ 
plant  the  parametric  continuity  constraints  when  build¬ 
ing  spline  curves,  the  bivariate  Beta  constraints  can  re¬ 
place  the  parametric  constraints  when  building  spline 
surfaces.  It  can  be  shown  that  the  Beta  constraint*  for 
first  and  second  order  are  equivalent  to  requiring  cob* 
tinuity  of  tangent  planes  and  Dupin  incat  rices  of  the 
patches  match  along  the  boundary  curve*.  Thus,  the 
chain  rule  approach  agrees  with  geometric  intuition  for 
both  Gl  and  C*  continuity.  Moreover,  the  chain  rule 
approach  yields  the  second  order  constraints  with  lees 
effort  than  the  geometric  approach.  For  higher  cider 
continuity,  geometric  intuition  becomes  more  feeble,  bat 
the  chain  rule  approach  still  applies. 

i  ■  ■  '  i 

4.  Conclusion 

We  have  defined  n,k  order  geometric  continuity  for 
parametric  curves  and  surfaces,  and  derived  the  Beta 
constraints  that  are  necessary  and  sufficient  for  it.  The 
derivation  of  the  Beta  constraints  is  based  on  a  simple 
principle  of  reparametrisation  in  conjunction  with  the 
univariate  chain  role  for  curves,  and  the  bivariate  chain 
rule  for  surface*.  This  approach  therefore  uncovers  the 
connection  between  geometric  continuity  for  carves  and 
geometric  continuity  for  surfaces,  provides  sew  insight  ' 
into  the  nature  of  geometric  continuity  in  general,  and 


'allows  the  determination  of  the  Beta  constraints  with 
lass  effor)  than  previously  required. 

The  use  of  the  Beta  constraints  for  G*  continuity 
allows  the  introduction  of  a  shape  parameters  for  curves, 
and  n(n  +  8)  shape  functions  for  surfaces.  The  shape 
parameters  and  shape  functions  may  be  used  to  modify 
the  shape  of  a  geometrically  continuous  curve  or  surface, 
respectively.  However,  geometric  continuity  is  only 
appropriate  for  applications  where  rat*  aspects  of  the 
parametrisstions  are  unimportant  since  discontinuities 

h  rate  are  allowed.  *  j  ,  -  • 

/ 

As  a  final  comment,  the  approach  we  have  taken 
is  not  based  on  measures  that  are  inherent  to  curves 
and  surfaces,  so  the  generalisation  to  4-variate  objects 
(volumes,  hyper-vohimes,  etc.)  can  be  made  very  sim¬ 
ply:  two  4-variate  parametrisatioas  are  GO- equivalent 
if  they  are  related  by  a  change  of  parametrisation  with 
positive  Jacobian;  the  corresponding  Beta  constraints 
may  be  derived  in  complete  analogy  to  the  development 
of  Section  S,  using  the  4-variate  chain  rule*  in  place  of 
the  bivariate  chain  rule. 
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